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ABSTRACT 

New observations from the Hubble ultra deep field suggest that the star formation rate at 
z > 7 drops off faster than previously thought. Using a newly determined star formation 
rate for the normal mode of Population II/I stars (PopII/I), including this new constraint, we 
compute the Thomson scattering optical depth and find a result that is marginally consistent 
with WMAP5 results. We also reconsider the role of Population III stars (PopIII) in light of 
cosmological and stellar evolution constraints. While this input may be needed for reioniza- 
tion, we show that it is essential in order to account for cosmic chemical evolution in the early 
Universe. We investigate the consequences of PopIII stars on the local metallicity distribution 
function of the Galactic halo (from the recent Hamburg/ESO survey of metal-poor stars) and 
on the evolution of abundances with metallicity (based on the ESO large program on very 
metal-poor stars), with special emphasis on carbon-enhanced metal-poor stars. The metallic- 
ity distribution function shape is well reproduced at low iron abundance ([Fe/H]> —4), in 
agreement with other studies. However, the Hamburg/ESO survey hints at a sharp decrease of 
the number of low-mass stars at very low iron abundance, which is not reproduced in models 
with only PopII/I stars. The presence of PopIII stars, of typical masses 30-40 M , helps us 
to reproduce this feature, leading to a prompt initial enrichment before the onset of PopII/I 
stars. The metallicity at which this cut-off occurs is sensitive to the lowest mass of the massive 
PopIII stars, which makes the metallicity distribution function a promising tool to constrain 
this population. Our most important results show that the nucleosynthetic yields of PopIII 
stars lead to abundance patterns in agreement with those observed in extremely metal-poor 
stars. This can be demonstrated by the transition discriminant (a criterion for low-mass star 
formation taking into account the cooling due to C II and O I). In this chemical approach to 
cosmic evolution, PopIII stars prove to be a compulsory ingredient, and extremely metal-poor 
stars are inevitably born at high redshift. 

Key words: nucleosynthesis, abundances - methods: statistical - stars: abundances, chemi- 
cally peculiar, formation - cosmology: miscellaneous. 



1 INTRODUCTION 

Recent data from the Hubble Ultra Deep Field (HUDF) have al- 
lowed the determinations of the cosmic star formation rate (SFR) 
to be extended from redshifts z ~ 4 up to z ~ 7 — 8 
jBouwens et al.ll2007h . In the context of the ACDM model, which 
has by and large been confirmed by the observation of anisotropies 
in the dif fuse cosmological microwave backgrou nd (CMB) by 
WMAP JSpergel et al] 120071 ; iDunklev et ail l2009h . a remaining 

* E-mail: rollinde@iap.fr (ER); vangioni@iap.fr (EV); dmau- 
rin@lpnhe.in2p3.fr (DM); olive@physics.umn.edu (KAO); daigne@iap.fr 
(FD); silk@astro.ox.ac.uk (JS); vincentf@iap.fr (FHV) 
f Institut Universitaire de France. 



challenge lies in understanding the first structures, including galax- 
ies and stars. This step requires either direct observations (of black 
holes, gamma-ray bursts. . . ) or indirect constraints through lumi- 
nosity functions , metal pollution and relic metal-poor stars (e.g. 
Tumlinson 2006; Salvadori et al. 2007). Our work develops the lat- 
ter approach, and focuses on the nucleosynthesis pollution from the 
first stars and their consequences for relic metal-poor stars. 

The role of massive PopIII stars as ionization sources at high 
redshift is poorly understood. Indeed, their role has been clouded 
by the fact that the measured value of the Thomson scattering op- 
tical depth (from z = to the redshift of emission of the CMB) 
has decreased significantly from WMAP1 to WMAP5. Current 
measurements imply a redshift for an instantaneous reionization 
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of z = 11.0 ± 1.4 (r = 0.087 ± 0.017) dPunklev et afll2009l) . 
In this context, the new observations of the cosmic SFR at high 
redshift bring a fresh perspective on the role played by PopIII 
stars. Since PopIII stars have a very specific impact on nucleosyn- 
thesis, it is useful to incorporate constraints from stellar obser- 
vations. From the nucleosynthetic point of view, halo stars have 
long been used to constrain galactic chemical evolution models, 
but were quite disconnected from cosmological models. Yet PopIII 
stars might be the primary source for the early metal enrichment 
of the interstellar medium (ISM) as well as for the intergalactic 
medium (IGM). The different abundance patterns observed in ex- 
treme meta l-poor stars (EMPS) may well be e xplained in terms of 
these stars dYoshida et alj|2007l ; lFrebelll2008j) . It has been shown 
( Brom m & Loebll2003h that the abundances of ionized carbon and 
neutral ato mic oxygen a re important for the transition from PopIII 
to PopII/I. IFrebel et all J2007I) have defined a transition discrimi- 
nant, Dtrans, and we show below that this quantity clearly reveals 
the nucleosynthetic imprint of PopIII stars. 

The paper is organized as follows. In Section [2] we summa- 
rize our model for the global chemical and cosmic evolution. New 
constraints related to the SFR and reionization are described in Sec- 
tion[3] We first reconsider the reionization constraint from WMAP5 
data in terms of a single normal mode of PopII/I star formation , 
by fitting the SFR to the most recent data l lBouwens et alj|2007h . 
We compare this result to a model which includes a contribution 
from PopIII stars. Stellar observations are described in Section|4] 
To address these issues, we use some of the abundances (Fe, O, C, 
Si) of individual stars of the Euro pean Southern Obse rvatory Large 
Program (ESO-LP) 'First Stars' JCavrel et af]|2004) and of pecu- 
liar extremely-low metallicity stars. We also compute the metallic- 
ity distribution function (MDF, number of observed low-mass stars 
at a specific metallici ty) and compare this to recent observations 
( Scho erck et al.| [2008). A grid study of the model parameters al- 
lows us to draw confidence contours for the SFR of PopIII stars 
and on the transition discriminant in Section [5] We also consider 
abundances derived from ID and 3D model atmospheres and com- 
pare results in each case. Our best model is used as an illustration 
in Sections[3]and|4] We summarize our results and conclude in Sec- 
tion!)] 



2 A GLOBAL COSMIC EVOLUTIONARY MODEL 

Any model for star formation in a cosmological context requires the 
inclusion of a model for dark matter structure formation, accretion 
and outflow of baryonic matter. Numerical simulations consistently 
follow the dark matter and baryonic components. Merger trees have 
also been used extensively. These models have the ability to probe 
different IMFs and critical metallicities. However, as yet, none of 
these methods follow individual element abundances. In a comple- 
mentary approach, we have developed a detaile d model of cosmo- 
logical chemical evolution (Daigne et al. 2006), using a simplified 
description of non linear st ructures, based on the s t andard Press- 
Schechter (PS) formalism JPress & Schechteri 1 19741: Ijenkins et ail 
l200ll ; lsheth &Tormenll 19991) . 

One of the main advantages over semi-analytical or fully re- 
solved simulations (that take up a lot of computer time) is the ability 
to probe a large region of parameter space and to follow abundances 
of individual elements (we use only Fe, C and O abundances in 
this study). Such simplified models are successfully used to tackle 
speci fic questions related to early star formation and reioniza- 
tion {Kampakoglou et al. 2008; Henri ques et alj [2009; Bagl a et all 



2009). These models only give average quantities, but they can re- 
produce all salient features of semi-analytical approaches. They are 
also supported by the lack of scatter in chemical observations of 
stars, and the fact that the MDF from Galactic halo field stars is 
statistically very similar to that of Galactic globular cluster sys- 
tems and th e stellar population of the nearest dSph satellites of 
the G alaxy iSchoerck et ai]|2008l ; iKirbv et alj|2008l ; IFrebel etafl 
2009). However, note that the Galactic globular cluster MDF drops 
to zero at [Fe/H]=-2.4. Presently, it is unclear where the cut-off is 
for the stellar population of dSph. 

We consider two distinct modes of star formation: a nor- 
mal mode of PopII/I and a massive mode of PopIII stars. The 
observed SFR and element abundances at redshift z < 6 (Sec- 
tion|4j are accomm odated by the normal mode with a Salpeter IMF 
toaigne et alj|2006h . while the WMAP5 integrated optical depth is 
marginally reproduced (Section [3}. For PopIII stars, we perform 
a grid analysis of a few key parameters to check their relevance 
for the chemical abundances of very metal-poor stars. Throughout 
this paper, a primordial power spectrum with a power law index 
n = 1 is assumed and we adopt t he cosmological parameters of the 
so-ca lled concordance model JSpergel et all 120071 : iDunklev et all 
l2009h . i.e. fi m = 0.27, Ov = 0.73, h = 0.71 and erg = 0.9. 

2.1 Description of the model 

Baryons are divided into three reservoirs. Two reservoirs account 
for the matter in collapsed structures: the gas (hereafter "ISM") and 
the stars and their remnants (hereafter "stars"). The third reservoir 
represents baryons in the medium between these structures (here- 
after "IGM"). The evolution of the mass of each reservoir is gov- 
erned by equations (1) and (2) in Daigne et al. 2006 (see also figure 
1 in the same paper), which depend on different fluxes: (i) baryon 
accretion from the IGM to the ISM, due to the structure formation 
process. This term is computed at each redshift in the framework 
of the Press & Schechter formalism (see equation (6) in Daigne 
et al. 2006), assuming that collapsed structures correspond to dark 
matter halos with a minimum mass M m i n (which, for simplicity, is 
kept constant over the entire evolution); (ii) baryon ejection from 
the ISM into the IGM, associated with outflows from the collapsed 
structures. Such outflows are powered by a fraction e of the kinetic 
energy released by SN explosions. The corresponding mass flux is 
evaluated by assuming that the typical velocity of these outflows 
scales as the mean escape velocity at redshift z. This velocity is 
computed by averaging the escape velocity over the distribution of 
the mass of dark matter halos above M m i n at redshift z. Outflows 
become less efficient with time as the typical mass of collapsed 
structures increases; (iii) the third flux is associated with the star 
formation process and transfers baryons from the ISM to stars. It 
can easily be computed once a star formation rate 3 , (z) has been 
specified. Our assumptions regarding this SFR are detailed below; 
(iv) the final flux stands for mass ejection from stars to the ISM, as- 
sociated with stellar winds and explosions. It is computed from the 
SFR and the IMF. We do not assume instantaneous recycling and 
therefore take into account the lifetimes of stars. The stellar data 
we use are described in the next subsection and are dependent on 
both the initial mass and metallicity of the stars. 

In addition to the evolution of the mass, we also compute the 
chemical evolution in the "ISM" and the "IGM". For any chemical 
element i, the mass fractions in the ISM and the IGM, Xj SM and 
X} GM , are followed as a function of the redshift. The correspond- 
ing differential equations can be found in Daigne et al. 2004 (equa- 
tions (6) and (7)) and are derived from the previous mass fluxes 
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assuming that (i) the chemical composition of the accretion flow 
related to the structure formation process is the composition of the 
IGM at redshift z; (ii) the chemical composition of the outflows 
escaping collapsed structures is the composition of the ISM at red- 
shift z; (iii) this composition is also the composition of stars that 
form at the same redshift; (iv) finally the composition of the matter 
ejected by stars into the ISM is given by stellar yields. These yields 
are detailed in the next subsection. Once the chemical evolution is 
known, it can be inverted to estimate the redshift of formation of a 
star from its observed abundance, as explained in Sec tion l5. 1.21 

The models considered here are bimodal. Each model contains 
a normal mode with stellar masses between 0.1 Mq and 100 Mq, 
with an IMF with a near Salpeter slope. The SFR of the normal 
mode peaks at z ~ 3. In addition, we allow for a massive compo- 
nent, which dominates star formation at high redshift. We primarily 
focus on stars with 40-100 M©, which terminate as type II super- 
novae. The nucleosynthetic pattern of pair-instability SN (PISN) is 
briefly commented on in the conclusions. 

In addition to the SFR, the IMF and stellar data, this simple 
model introduces only two other parameters: the minimum mass 
of the dark matter halo in collapsed structures, M m i n , and the ef- 
ficiency of supernovae to power outflows from collapsed struc- 
tures, e. The impact of these two parameters has been discussed 
in detail in Daigne et al. 2006. In the present paper, we adopt 
a typical scenario with M m in = 10 7 Mq and e = 3 x 10 -3 . 
This minimum mass, M m i n , is also d educed from h y drody nami- 
cal simulations at z ~ 15. Note that Ijohnson et al.l (120091) con- 
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sider M mi n = 10 8 M B at z = 12.5 (see also iDaigne et ai]|2006l : 



iBagla et alj|2009l ; iTrenti & Stiavellill2009h . e primarly impacts on 
the enrichment of IGM, which is not considered in this paper. 



2.2 Yields and lifetimes 

The lifetim es of intermediate mass s tars (0.9 < M /Mq < 8) are 
taken from lMaeder & Mevnetl dl989l) and from lSchaereiH2002h for 
more massive stars. Old halo stars with masses below ~ 0.9 M© 
have a lifetime long enough to be observed today. They inherit the 
abundances of the ISM at the time of their formation. Thus, their 
observed abundances reflect, in a complex way (due to exchanges 
with the IGM), the yields of all massive stars that have exploded 
earlier. 

The yields of stars depend on their mass and their metallic- 
ity, but not on their status (i.e. PopII/I or PopIII). Some PopII/I 
stars are massive, although in only a very small proportion since 
we use a steep Salpeter IMF. PopIII stars are all very massive 
sta rs. We use the tables of yields (an d remnant types) given 
in Ivan den Hoek & Groenewegenl d 19971) for intermediate m ass 
stars (< 8Mq), and the tables in IWooslev & WeaveJ J 19951) for 
massive stars (8 < M/Mq < 40 ). An interpolation is made 
between different metallicities (Z=0 and Z=10" 4 ' ~ 3 ' ~ 2 ' _1 '°Z©) 
and w e extrapolate the t abulat ed values beyond 40 Mq . 

Heger & Wooslevl d2008l) have provided new stellar yields in 
the mass range 10 — 100 M©, but only at zero-metallicity. At zero- 
metallicity, and in the 40 — 100 Momass range, we have checked 
that our extrapolated values are consistent with those yields. In- 
deed, their favored model corresponds to very little mixing and a 
high o xygen to iron ratio, in agreement with the lWooslev & Weaver! 
(1995) model we chose. We expl icitly checked that o ur res ults 
were not modified usin g either IWooslev & Weaver! dl995l) or 
iHeger & Wooslevl d2008l) at zero-metallicity. 



0.1 



m 10 ■> 



0.5 
I 



0.1024 
L 



-Pop [I/] 




Figure 1. Cosmic SFR as a function of redshift. PopII/ I: The observed SFR 
up to z ~ 5 (solid blue points) is taken fromlHop kins & Beacoml J2006h . 
The dashed blue data point comes from Bouwens et al. (2007). The solid 
black line is the extrapolated PopII/I SFR, using the shape given by Eq. (3, 
in order to match the data at z ^ 5 and the upper limi t at z ~ 7. PopIII : 
The recent non-observation of Lyo-Hell dual emitters jNagao et al . 2008) 
at z ~ 4 — 4.5 is shown with the red arrow. The red line corresponds to 
the PopIII SFR upper limit, adjusted to satisfy the Icr WMAP5 upper limit 
on the Thomson optical depth (see Fig. |2). This limit has to drop (dotted 
line) to be reconciled with the Lyo-Hell observation. The envelope (shaded 
area) contains all models, that reproduce the Dtrans evolution at a 95% CL 
(see Section \5\ and Fig. |4) For illustration, the best model for the PopIII 
SFR (including nucleosynthetic and cosmological constraints, described in 
Section |4) is shown by the dashed red line. This SFR is parametrized with 
Eq. {T}, with the following parameters: um = 0.0016 Mq yr _1 Mpc~ a 
(astration rate for PopIII stars), z m ni = 22.8, am = 4 and bni = 3.3. 



3 COSMOLOGICAL CONSTRAINTS 

Takin g into account the res ults of recent observations on the PopII/I 
SFR dBouwens et alJ l2007). we reanalyse the reionization capabil- 
ity of these stars. The WMAP5 upper limit on the Thomson optical 
depth is then used to set constraints on the PopIII SFR. 



3.1 PopII/I SFR 

We fit the SFR his t ory of PopII/I stars to the data compiled in 
Hopki ns"& Beacoml d2006l) ( from z = to 5), and to the recent 
measurement at z ~ 7 by iBouwens et alj d2007l) (see Fig. [TJ . 
The behaviour of the SFR at high redshift is still disputed. Some 
indirect measurements hint against a decline of the SFR beyond 
z ~ 3 dFaucher-Giguere et al. 2008), but direct observations of 
■Z850 -dropou t galaxies seem to confirm a sharp decr ease of the SFR 
at z ~ 6—7 dOesch et alj2009h . The low level of the lBouwens et al.l 
d2007l) data thus places strong constraints on the PopII/I SFR. 

The fit for the SFR, ijj(z), is based on the mathematical form 
proposed in ISpringel & Hemquistl ( l2003l) : 



ip(z) = v- 



aexp(6 (z 



.)) 



(1) 



a — b + &cxp(a (z — z m )) 

The amplitude (astration rate) and the redshift of the SFR max- 
imum are given by v and z m respectively, while b and b — a 
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are related to its slope at low and high redshifts respectively. In 
the following, we use the subscripts II/I and III for parameters 
related to PopII/I and PopIII SFRs respectively. The thick black 
curve in Fig.Q]is computed with vxxn = 0.3 Mq yr _1 Mpc -3 , 
^mii/i = 2.6, dir/i = 1.9 an d bn/i = 1.2. The para meter values 
chosen here differ from that in iGreif & Bromml J2006h so as to bet- 
ter account for the high redshift data point dBouwens et alj|2007l) . 
To be conservative with respect to the reionization constraint (see 
below), we choose a SFR consistent with the la upper limit of this 
z = 7 observation. A steeper slope only decreases the number of 
ionizing photons. 



3.2 Reionization from PopII/I stars 

Having set the SFR in our model, we now compute the electron 
scattering optical depth for our choice of IMF. The evolution of the 
volume filling fraction of ionized regions is given by: 

dQiontz) 1 dn io n(z) 



dz 



dz 



ctBnbC(z) 



xQ? on (z)(l + z) 3 



(2) 



where rib ls the comoving density in baryons, ni on (z) the co- 
moving density of ionizing photons, «b the recombination coef- 
ficient, and C(z) the clumping factor. This factor is taken from 
IGreif & Bromml ( 2006) and varies from a value of 2 at z 20 to a 
constant value of 10 for z < 6. The escap e fraction, / esc , is set to 
0.2 for both PopIII and PopII/I SN, while lGreif& Bromml 120061) 
take two different values. The number of ionizin g photons for m as- 
sive stars is calculated using the tables given in Schaerer (2002). Fi- 
nally, the Thomson optical depth is computed as in IGreif & Bromml 
(2006): 



f 

Jo 



dz' Qion(z') (1 + z')" 



dt 



dz 1 



(3) 



where z is the redshift of emission, and ctt the Thomson scattering 
cross-section. 

In the upper panel of Fig. [2] we show the predicted evolu- 
tion of the volume filling factor of ionized regions, Qi on as a func- 
tion of redshift. In the lower panel, we show the integrated optical 
depth from z = to z . These results can be compared to Fig. 5 of 
IGreif & Bromml d2006h . We see that the normal mode (black solid 
curve) is contained in the la WMAP limit (grey band) for the ob- 
served value for the Thomson optical depth. 



3.3 Reionization from PopIII stars 

A massive PopIII mode at high redshift makes it easier to reproduce 
the optical depth at high redshift. For illustrative purposes, the his- 
tory of reionization is shown in Fig.|2]for a typical SFR (red dashed 
line in Fig.[TJ found to best meet the local constraints (discussed in 
Section[5j- The parameters of the best fit model for PopIII SFR are 
given in the caption of Fig. Q] The associated history of reioniza- 
ti on is consistent with a two-step model of reionization as invoked 
in Dun klev et all §009), who also comment that the WMAP5 data 
suggests a more gradual process with reionization beginning per- 
haps as early as 20. In particular, the PopIII SFR allows Qi on 7^ 
at z > 10. The bimodal SFR is consistent with recent general sce- 
narios that have emerge d from hydrodynamic simulations and ana - 
lytical calculations (e.g. lBromm & LarsorfcOOllGreif et al]. 2008). 
These studies have shown that, due to the re-collapse of the struc- 
tures after disruption by massive SNe, low-mass star formation may 
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Figure 2. Histories of volume filling fraction of ionized regions, Qi on (top 
panel) and the implied optical depth (lower panel) with and without PopIII 
stars (dashed red and solid black respectively). The envelope in the lower 
panel (light shaded area) shows the maximum and minium values of optical 
depth reached by the different models that reproduce the chemical abun- 
dance constraints at a 95% CL (see Section|5J. The lower panel also shows 
the range of a llowed values for the optical depth as measured by WMAP5 
iDunklev et alj2009l) . 



have been delayed by more than a Hubble time or up to a red- 
shift ~ 10. This can be understood as a statistical d escription of the 
mech anical feedback implemented in merger trees dSalvadori et al.l 
l2007h where a time delay between PopIII and PopII/I appears at 
each individual branch. 

Using a grid analysis on the parameters of Eq. l[T) for PopIII 
stars (see details in Section|5j, we are able to derive an upper limit 
on the SFR. There is no lower limit as PopII/I stars alone are con- 
sistent with the WMAP reionization lower limit. The 68% CL up- 
per limit on the SFR is reported as a solid red line in Fig. Q] that 
we stop below z < 10. At lower redshift, the most stringent con- 
strai nt is set by the rece nt non-observation of Lya-Hell dual emit- 
ters dNagao et al.ll2008l) , which gives an upper limit for the PopIII 
SFR of 5 • 10" 6 M yr _1 Mpc -3 at z ~ 4 - 4.5 (red arrow). 
The dotted line reconciles this observational constraint with the 
WMAP5 upper limit. Note that within the region allowed by the 
cosmological constraint (bottom right), the shape of the SFR of 
PopIII remains completely unspecified (based on the reionisation 
constraint), as the relevant quantity is related to the integral over 
z of the SFR. We show in the following section that local obser- 
vations require the existence of PopIII stars, and set stronger con- 
straints on the SFR and its lower limit. 



4 LOCAL OBSERVATIONS 

The average mass of the structures in the PS formalism is close to 
the total mass of the Milky Way (10 12 Mq) at z = 3, which corre- 
sponds to the epoch of the peak of the global SFR. Our average de- 
scription of the cosmic history is thus expected to be representative 
of the evolution of the Galaxy. In this section, local observations 
are analyzed in a global cosmological context. 
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Figure 3. Comparison of the observed MDF iSchoerck et al. 2008, thick 
histogram) for the galactic halo stars with model calculations of the normal 
mode of star formation (black solid line). Other lines show the effect of the 
inclusion of PopIII stars with different mass ranges. All distributions are 
normalized to match the observations at [Fe/H]=-3. 



4.1 Metallicity Distribution Function (MDF) 

The MDF is the distribution of stars at a given iron abundance, 
derived from observations that have been corrected for various 
technical and observationnal biases. It has been recentl y pub- 
lish ed by the Hamburg-ES O survey dSchoerck et al.ll2008h using 
the IChristlieb et al.| |2008l) catalog. The thick line histogram in 
Fig.l3lis taken from lSchoerck et aljd2008l) . As noted there, in order 
to compare observed and predicted MDF, one has to take into ac- 
count the modification of the shape of the MDF by the selection of 
metal-poor candidates. Consequently, a comparison to our predic- 
tion is valid for [Fe/H]< -2.5 onljQ. 

In our analysis, for a given model (SFR and IMF of PopII/I 
and PopIII), we count the number of low mass stars created at each 
redshift. As the iron abundance in the ISM is also calculated as a 
function of the redshift, the metallicity of the stellar population at 
each redshift is known (Section [2.1t . The value of the MDF in a 
given iron abundance bin, is given by the number of stars that are 
still shining today and were created at times when the iron abun- 
dance was within this bin. The MDF for several possible PopIII 
star formation histories are plotted along with the observation in 
Fig. [3] The solid black line corresponds to the normal mode, i.e. 
a Salpeter IMF starting at 0.1 Mq. A massive mode is then added 
with a Salpeter IMF with different lower masses: 36.5 Mq (red 
dashed), 30 M (black dotted) and 8 M (black dot-dashed). The 
minimum mass of 36.5 corresponds to our best model. Note that 
the value of the upper limit of the IMFs is not a key parameter, in 
constrast to the lower one, due to the steep slope of the IMF. 

For -4 <[Fe/H]< -2.5, the slope of the MDF is well 
reproduced by PopII alone. The reality of the drop observed at 
[Fe/H]~ —3.6 (compare the thin solid black line to the thick black 



histogram) is still debated (e.g., iKomiva et allbOOg^ The predic- 
tion of several models for the MDF is discussed in lschoerck et al] 
(2008), where the authors conclude that none of them are able to ex- 
plain the_tailat [Fe/H]< —4, except a stochastic enrichment model 
(e.g. lKarlssonl l2006). In our model, as in Karlsson's model, the 
presence of a massive mode implies the early production of iron, 
which has an impact on the low metallicity MDF: very few low 
mass stars at zero or quasi-null metallicity are expected in such a 
scenario, i.e., effectively as a prompt initial enrichment (PIE). The 
cut-off on the MDF at low metallicity depends on the lower mass 
of the massive mode. 

Given the paucity of stars with [Fe/H]< —3.5, it is difficult 
to accurately discriminate between different star formation his- 
tories presently. Indeed, PopIII stars do not modify the MDF at 
—3.5 <[Fe/H]< —2.5, as long as their minimal mass is larger than 
about 30 M Q , 

4.2 Global nucleosynthetic evolution 

We now focus on the giants of the ESO-LP dCavreletal .1120041) 
which contain the abundances of 17 elements from C to Zn for 
all of the observed s t ars. F ollowing the nomenclature proposed 
in iBeers & Christliebl d2005l) . very metal-poor (VMP) stars cor- 
respond to a metallicity [Fc/H] < —2, extremely metal-poor 
(EMP) stars to [Fe/H] < -3, ultra metal-poor (UMP) stars to 
[Fe/H]< —4 and hyper metal-poor (HMP) stars to [Fe/H] < —5. 
Some UMP stars are also carbon-enhanced metal-poor (CEMP) 
stars with [C/Fe]> +1. 

4.2.1 ESO-LP and CEMP data 

The ESO-LP contains 35 VMP stars, which includes 22 EMP stars. 
The star in their sample with the next to lowest iron abundance 
(CS 229 49-037; [Fe/H] = -3.97 ) is a CEMP star. Note that CD- 
38°245 teessell&NorrisIl 19841) has [Fe/H]=-4 19; unfortunately, 
ICavrel et al] d2004l) and lBessell & Norrisl d 19841) do not provide O 
abundances, and only an upper limit to the abundance of C. One 
result of this survey has been the reduction of the scatter in different 
element abundances as correlated against [ Fe/H] (e.g Mg, Ca, Cr 
and Ni). According to IBeers & Christliebl {2005), this is contrary 
to the long-standing hypothesis that, at such low metallicity, one 
observ es the nucleosynthetic products of only a few or even a single 
SN II jShigevama &Tsuiimotdl 19981 ; iTsuiimoto et al.ll2000h . The 
lack of scatter in these abundances could be explained by a well- 
mixed ISM. 

The CE MP stars re f erence s are retrieved from the SAGA tooQ 
presented in lSuda et"a 1. (2008). Abundances derived from ID and 
3D model atmospheres and references are gathered in Table Q] At 
the lowest metallicity, HE 0107-5240 and HE 1327-2326 are both 
HMP and CEMP stars. The very specific chemical pattern of these 
stars is assumed to be related to the first stages of star formation. 
The UMP star HE 1300+0157 with [Fe/H]=-3.88 has an abundance 
pattern close to the single UMP star in the ESO-LP. For G77-61, the 
criterion of CEMP-no star is not strictly met (only an upper limit 
exists on [Ba/Fe]). Moreover, in contrast to the other stars consid- 
ered, the abundances of G77-61 are based on medium resolution 
optical ([Fe/H], [C/Fe]) and near-infrared ([O/Fe]) spectroscopy. 
This method can br ing systematic un certainties and possibly an 
over-estimate of C dBeers et all 120071) . Consequently, this star is 



1 For a given element X, the abundance relative to the solar value is defined 
by[X/H] = log 10 (X/H) - log lo ((X/H) ). 



http://saga.sci.hokudai.ac.jp/Retrieval/db.cgi 
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Table 1. Abundances of CEMP stars 

Name [Fe/H] [C/H] [O/H] 

1D/3D* 1D/3D 1D/3D 



HE 1300+0157 a 


-3.88/- 


-3.88 


-2.50/- 


-3.10 


-2.12/- 


-2.72 


HE 1327-2326 6 


-5.46/- 


-5.96 


-1.55/- 


-2.18 


-1.76/- 


-2.54 


HE 0107-5240 c 


-5.30/- 


-5.50 


-1.30/- 


-2.30 


-3.00/- 


-3.60 


G-77-61 d 


-4.00/- 


-4.20 


-1.40/- 


-2.40 


-2.20/- 


-2.80 



Abundances derived from ID and 3D model atmospheres * Abundances 
derived from 3D/NLTE model atmospheres are used for the main analysis 
(Section [5}. Abundances derived from ID model atmospheres are consid- 
ered in Section l5.3.2l 

" assumed to be a s ubgiant as in iFrebel et aljj2007l) . 
b |Frebel et al.H2008l) 

c lBessell et al.H2004l) ; lchristheb et alj j2004l); abundances derived from 3D 
model atmospheres are estimated from Coll et et al . (2006). 
d iBeers et al] l2007l) : ICollet et alj {iwSi . 

not used in the analysis, but is displayed in the figures for illustra- 
tive purposes. There is one additio nal star with —5 < [Fe/H] < —4, 
HE 0557-4840 jNorris et al.l2007t) . While it also has an enhanced C 
to Fe ratio, there is only an upper limit on [O/H] making it difficult 
to place this star directly into our analysis. 

For our main analysis, abundances derived from 3D model at- 
mospheres are considered first. We also discuss the impact of using 
ID versus 3D model atmospheres in Section [5.3.2l 

4.2.2 Transition discriminant 

The transition from one star formation mode to another is of- 
ten described in term s of a critical metallicity. For example, 
iBromm & Loebl J2003h have shown the importance of ionized car- 
bon and neutral atomic oxygen in this transition and suggest that 
when sufficiently abu ndant, these eleme nts act as a trigger to lower 
mass star formation. IFrebel et alj J2007t) defined the transition dis- 
criminant: 

Dtrans = log 10 (10 [C/H] + 0.3 X 10 [ ° /H] ) . (4) 

It is argued that low-mass star formation require s Dtrans ~ — 3.5 ± 
0.2 dSantoro & Shullll2006l : IFrebel et al.ll2007ri . A scatter plot of 
Strans as a function of [Fe/H] is shown in Fig. [4] All ESO-LP 
very metal-poor stars (circles) and CEMP stars meet the condition 

Dtrans ~ 3.5. 

4.2.3 PopIII stars and the enrichment of CEMP stars 

Our model follows C, O and Fe abundance evolution with z, and 
hence enables us to calculate -Dtrana and [Fe/H] at each redshift 
(see Section [2] for details and iDaigne et alj|2006t) . In the -Dtrans 
versus [Fe/H] scatter plot of Fig. [4] the evolution from z — 30 
to 2 = corresponds to a continuous path starting at the lowest 
metallicity (left-hand side) and ending up at the highest metallicity 
(right-hand side). 

A standard IMF used for the normal mode leads to the thick 
black-solid line. It reproduces well the ESO-LP observations, but 
it is clear that the CEMP stars abundances are not explained in this 
standard scenario. The evolution of Dtrans in our best fit model 
(see next section) which includes a massive mode is shown as the 
red-dashed line. The SFR corresponding to this model is shown in 
Fig. Q] (red-dashed line) and corresponds to the parameters given 
above in Section [3~3l 
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Figure 4. The transition discriminant, Dtrans, for metal-poor stars as a 
function of [Fe/H]. The ESO-LP stellar data are shown as circles while pe- 
culiar CEMP stars are indicated with different symbols. The thick black 
solid line corresponds to a model with PopII/I stars alone. The addition 
of PopIII stars (35-100 Mq, our best model) leads to the red dashed line. 
The dotted vertical line indicates the prompt initial enrichment (PIE) pro- 
duced by PopIII stars (see Fig. [5J. The envelope (thin black curves) con- 
tains all evolution curves predicted by models, that match the abundance 
at a 95% CL, as determined by our x 2 analysis (see Sectionf5). Only stars 
with error bars are included in this analysis. 

PopIII SN eject more C and O than PopII/I yielding higher val- 
ues of Dtrans - This allows one to explain the abundances of the pe- 
culiar stars. As a result, these stars must be associated with star for- 
mation occurring within regions havi ng a specific nu cleosynthetic 
history at early times as suggested in iFrebell |2008a). At interme- 
diate redshift, we see a decrease of both -Dtrans and [Fe/H]. This 
corresponds to the epoch where the SFR is reduced (see Fig. [TJ. 
The accretion of pristine IGM material continues, leading to the di- 
lution of ISM abundances. The abundance dilution allows one to 
reproduce simultaneously the observations of CEMP stars at high 
redshift and of ESO-LP stars at low redshift. As PopII/I stars turn 
on, we reconnect to the black solid line. Finally, the vertical dotted 
line indicates the PIE produced by PopIII stars as discussed already 
for the MDF drop in Section HTI 



5 CONSTRAINTS ON POPIII STARS FROM LOCAL 
OBSERVATIONS 

As shown above, PopII/I stars barely account for the reionization 
of the Universe. PopIII stars are not essential in this context, un- 
less the SFR of PopII/I stars at z > 7 drops off faster than in our 
conservative choice (see Fig.[T|l, or the value of the escape fraction 
used in our study (20%) is overestimated. Conversely, if we turn to 
the local observations and interpret them in a cosmological context, 
PopIII stars are mandatory from chemical observations, especially 
based on the consideration of CEMP stars. Hence, if we consider 
that the abundances of CEMP stars are representative of the nucle- 
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osynthesis of the first stars, we can perform a statistical analysis to 
constrain the PopIII parameters. We have shown above that the few 
stars observed at very low metallicities are not sufficient to use the 
MDF as a constraint on PopIII models. To reproduce the slope of 
the MDF at [Fe/H]> —3.5, we consider only minimal masses of 
30-40 Mq for PopIII stars. 

We first describe the method to evaluate the quality of the 
model. In this process, a redshift of formation can be set to each 
observed star. We then discuss the constraints set by abundances 
derived from 3D model atmospheres and the uncertainties related 
to ID or 3D model atmospheres. 

5.1 Methodology 

The analysis is done in the D t rans-[Fe/H] plane, using a grid ap- 
proach to cover the parameter space of the model. The parameters 
that we vary are Uin (from 2 x 10" 5 to 2 x 10~ 2 M yr" 1 Mpc -3 ), 
z m iii (from 10 to 26), am (from 0.5 to 5), 6 m (from 0.03 to 4.9) 
and the minimum mass of the massive mode, Mm (from 30 to 39 
Mo). 

5.1.1 Minimization procedure 

For each model, we calculate the D tra ns-[Fe/H] curve, which is 
multi- valued for models including PopIII stars. The best fit model 
corresponds to the set of parameters that minimizes a x 2 on the 
predicted and observed stellar iron abundance and transition dis- 
criminant Dtrans (via stellar carbon and oxygen abundances) of a 
sample of 21 ESO-LP and 3 peculiar stars. We consider all param- 
eters within a 95% confidence level (CL) based on this \ 2 to be 
acceptable. Therefore, in all figures, we have displayed 95% CL 
envelopes that contain the predicted curves of all those acceptable 
models. 



5.1.2 Transforming stellar abundances to redshift 

In many models, the iron abundance is monotonically increasing 
with time, giving a unique correspondance between [Fe/H] and red- 
shift. This is the classical view of cosmological chemical evolution, 
where the oldest stars are the most metal-poor. However, in the first 
metal-free structures, the nucleosynthesis of PopIII stars provide 
features in the evolution of [Fe/H], so that the previous picture does 
not necessarily hold. In that case, not only the abundances of Fe, but 
also of elements such as C and O are required in order to be able to 
assign the VMP stars to a particular redshift. 

Each position along the predicted curve D t ran S /[Fe/H] is as- 
sociated with a redshift. Following the procedure described above, 
we search for the position along this curve that minimizes the dis- 
tance to the observed values for a given star. This position gives the 
redshift of formation of this star. The predicted curve depends on 
the model considered, and therefore different models will associate 
different redshifts of formation with the same star. As we repeat this 
procedure for models that satisfy the 95% CL on the whole sample, 
we can infer the allowed range of redshift (from the minimum to 
the maximum redshift) for one specific star. 

5.2 abundances derived from 3D model atmospheres 

Models that satisfy a 95% CL limit are included within the enve- 
lope shown in Fig. [4] (thin black lines). All stars included in our 
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Figure 5. Evolution of the [Fe/H], [C/H] and [O/H] abundances with red- 
shift. Dashed red curves and black solid lines correspond to the prediction 
from the best model and to the 95% CL envelopes (SectionfS). The SFR of 
the best model is shown in Fig.[T] Symbols used for the stars are the same 
as in Fig. [4] Horizontal lines in the bottom panel give the estimated errors 
on redshift for CEMP stars (Tableff). 

analysis are located along the path defined by the best fit. As ex- 
plained in Section l4~2l G77-61 displays larger errors on its abun- 
dance determination. Yet, it is also located within the envelope 
(square in Fig.©. HE 0557-4840 is located at a somehow low value 
for [Fe/H]. According to the exact value for Dtrans, it could either 
seat outside the envelope (Dtrans— -2.2) or be very close to the edge 

(^trana— -3). 

5.2. 1 Consequence for SFR and optical depth 

The envelope of the PopIII SFR for all allowed models is shown 
in Fig. Q] (shaded area). Some common features can be noted: the 
massive mode becomes sub-dominant for redshift lower than about 
15; on the contrary, it is required at a level of 2x 10 -4 - 3 x 10 -3 
Mq yr -1 Mpc" 3 around 2 ~ 24 - 28. The PopIII SFR cannot 
be larger than 3 x 10~ 3 Mq yr -1 Mpc -3 independently of the 
redshift. 

The envelope of the optical depth history among all models is 
shown in Fig. [2] (light shaded area). We verify that all models are 
compatible with the la limit of WMAP5. 

5.2.2 Redshift of formation 

As explained in Section [5.1.2| for a given model, we have assigned 
a typical redshift, 2*, of formation to each individual star. Thus, 
we can locate each star on a plot of abundances versus redshift, 
and compare them directly to the predicted evolution. This is done 
in Fig. fusing the best fit model (dashed red line). We also show 
the corresponding 95% CL envelopes for the evolution of the abun- 
dances (thin black lines). Those envelopes also bracket the errors 
on the redshift of formation (see [Fe/H] plot in Fig. [5}. The redshift 
ranges are given in Table|2] 

We see that CEMP stars are likely to form at high redshift. For 
example, the star HE 1300+0157 (upper triangle in Fig.[4} is found 
to be formed at 2* ~ 11 (for the best fit model). HE 1300+0157 
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is located along a branch of the evolution where the dilution by 
IGM accretion is the only important process. Then, abundances de- 
crease very slowly and uncertainties on the redshift are of the order 
of 100 Myr (15 < z < 20). It is interesting to note that the pre- 
dicted abundances also reproduce the Fe, C and O abundances at 
this specific redshift. The star G 77-61 (square), that has a similar 
iron abundance but a higher carbon abundance, is born earlier (it is 
located on a different place along the curve .Dtrans -[Fe/H]; Fig.|4j 
at z ~ 16. The exact redshift for HE 1327-2326 and HE 0107-5240 
(lower triangle and diamond) depends on the exact time of forma- 
tion of the first structure, which is largely unknown. The duration 
of this episode is only a few tens of Myr, which also corresponds to 
the uncertainty on the redshift of formation of those stars. Finally, 
all stars observed in the ESO-LP are found to form at redshift be- 
low 10, except CS 22949-037 whose abundances are very similar 
to HE 1300+0157 (Sectionl4~2l. 



5.3 Uncertainties 

5.3.1 Peculiar stars 

The observed oxygen abundance in the UMP star HE 0107-5240 
(diamond in upper panel of Fig. [5} seems to be lower than the 
predicted abundances. This discrepancy points to the very special 
stars where C, O abundances are not explained b y standard predic- 
tions dNomoto et alj2006! ; lTominaga et al.l2007l) . It is clear that the 
yields of the massive stars at zero metallicity are poorly known. For 
example, new calculations concerning the effect of rotation on the 
evolution of primordial stars can modify consider ably the abun- 
dance of C and O, and their corresponding ratio fekstrom et al.l 
2008). They would then allow for a decrease of about 1 dex re- 
quired for the oxygen abundance of HE 0107-5240 to be repro- 
duced. 



5.3.2 ID vs 3D model atmospheres 

We now study the influence of uncertainties on stellar abundances 
on the overall results. We have performed the same full analysis 
using the abundances derived from ID model atmospheres (given 
in Table [TJ, yielding new envelopes for the _D t rans-[Fe/H] plane, 
abundance evolution, optical depth evolution and SFR of PopIII 
stars. In the following figures, the same symbols as above are used 
for each CEMP star. We do not show ESO-LP stars in the figures, 
as their positions are unchanged. 

The main effect of the correction applied to abundances de- 
rived from ID model atmospheres is to lower all abundances by 
about 0.1 dex for [Fe/H], and up to 1 dex for C and O abundances. 
This can be seen in Fig. [6] by comparing filled and open symbols 
(abundances derived from 3D and ID model atmospheres respec- 
tively). Given the predicted path of -Dtrans -[Fe/H], this does not 
modify much the envelope in the intermediate redshift range. At 
high redshift, the effect is strong and displaces the envelope towards 
lower values of Dtrans. In both cases, PopIII stars are required and 
a very similar pattern is displayed. 

As for the SFR, the best fit model assuming abundances de- 
rived from ID model atmosphere extends slowly down to about 
10~ r Mq yr _1 Mpc~ 3 at z ~ 4. However, the global envelope 
(including all possible shapes) is almost identical in both cases, al- 
though a bit more extended assuming ID model atmospheres. The 
evolution of optical depth is also not modified in a significant way, 
since the contribution of PopIII stars to reionization is marginal. 
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Figure 6. The transition discriminant, Dtrans, as a function of [Fe/H] (same 
as Fig.|4j. The abundances of CEMP stars derived from ID and 3D model 
atmospheres (Table [T} are marked with filled and open symbols respec- 
tively. The 95% CL envelopes of predicted evolution assuming 3D and ID 
model atmospheres are displayed with thick and thin solid lines respec- 
tively, while best models are shown with thick and thin dashed lines. 



Table 2. Estimation of the redshift of formation of CEMP stars 

Name 3D model 6 ID model 6 

Z n\in ^bcst 2 max ^min ^bcst ^max 



HE 1300+0157 


10.4 


11.6 


11.9 


9.7 


11.5 


12.0 


HE 1327-2326 


26.1 


27.9 


28.3 


24.9 


28.3 


28.4 


HE 0107-5240 


25.4 


27.6 


28.2 


24.7 


28.3 


28.3 


G-77-61 


14.3 


15.8 


18.0 


14.3 


15.9 


19.5 



Minimum and maximum redshift predicted by models that match 
abundances of all stars at a 95% CL (Section|5J. b Abundances are derived 
from 3D model atmospheres (second to fourth columns) or from ID model 
atmospheres (fifth to last columns). 



Finally, we consider how our estimates for the redshift of for- 
mation of CEMP stars are affected. Fig.|7]displays the envelope for 
the evolution of abundances with redshift allowed in the two cases. 
It appears even more clearly that the discrepancy occurs only at 
high redshift, while the transition from CEMP to standard metal- 
poor stars is required in both cases. The stars are not located in this 
figure for clarity. The different ranges of estimates for their redshift 
of formation are summarised in Table|2] Again, only the two stars at 
the highest redshift are affected by this uncertainty; while the range 
of allowed redshift is larger in the case of abundances derived from 
ID model atmospheres (this is not visible in Fig.|7)- 



6 DISCUSSION AND CONCLUSIONS 

Following the approach developed in iDaigne et alj fc004h . we have 
modeled the evolution of individual element abundances in the ISM 
assuming homogeneous star formation and s tellar yields. Recent 
observations at z~ 7-8 teouwens et alj[2007h were used to better 
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Redshift 

Figure 7. Evolution of the [Fe/H], [C/H] and [O/H] abundances with red- 
shift (same as in Fig. [5}. The 95% CL envelopes of allowed evolutions as- 
suming abundances derived from 3D and ID model atmospheres for CEMP 
stars are displayed with thick and thin lines respectively. The envelope as 
well as the evolution for the best models (thick and thin dashed red lines) 
are quite different at high redshift only. 



constrain one ingredient of the model, namely the SFR for PopII/I 
at high redshift. We have shown that a homogeneous scenario of 
hierarchical structure formation reproduces many different obser- 
vations, from reionisation and first star abundances to local abun- 
dance observations. We found that using the most recent results 
on the optical depth from WMAP, a massive mode is not abso- 
lutely required. Nevertheless, the data can accommodate a PopIII 
contribution, responsible for the gradual reionization starting from 
z ~ 20. Although the cosmological importance of PopIII stars can- 
not be fully constrained by the integrated Thomson optical depth, 
this question may be better tackled in the future with the help of ac- 
curate measurements of the CMB polarization data ( Dunkl ev et al.l 
120091) . 

We have also considered stellar constraints, in particular the 
MDF and the evolution of Dtrans with [Fe/H]. In the literature, the 
slope of the MDF can be reproduced by different models of chemi- 
cal enrichment, from galactic models to hierarchic al analytic mod- 
els inc luding me rger trees, as show n in Fig. 12 of ISchoerck et al.l 
(2008) (see also lTmnimsorJ l2006). This indicates that the slope 
of the MDF does not discriminate between the methods used and 
the associated level of heterogeneity. In c ontrast, the pattern a t 
[Fe/H] < —3.5 is more difficult to reproduce. ISchoerck et al.l ( 2008) 
claim that no model considered in their paper can reproduce the 
MDF at very low iron abundances. We find that the modification 
of the tail at low values of [Fe/H] may be related to the presence 
of PopIII stars. In our study, we have tested the impact of the mas- 
sive mode by varying the typical PopIII minimal mass from 8 to 
40 M© . This appears to be an important parameter as found in the 
study of Salvado r! et alj J2007I) . They used a mass range from 140 
to 200 M© which we consider disfavoure d by the very spe cific 
yields of the stars within this mass range toaigne et alj|200q) . In 
this paper, we demonstrate that the minimum mass must be larger 
than about 30 M© in order to reproduce the observed part of the 
MDF. 

We have shown that the existence of PopIII stars at high red- 



shift is required to explain the abundance pattern observed in the 
CEMP stars. In addition, we have shown that a massive mode with 
a typical mass of 40 Mg reproduces the evolution of observed 
D trans- In contrast, it is known that PISN with masses 140-200 Mq 
do not provide the correct chemical pattern and cannot reach high 
values of -Dtrans for low values of [Fe/H] (the ratio C/Fe in their 
yields, iHeger et alfcool is not high enough). However, if it could 
be established that th e CEMP stars were particular case s, such as 
belonging to binaries (Ryan et al. 200?]; lTumlinso n 2007) or due t o 
the preferential depletion of iron in grains JVenn & Lambertl2008l) . 
the SFR related to PopIII stars would be diminished, at least as far 
as chemical evolution is concerned. 

In conclusion, our analysis hints at a massive mode at 
z ~ 20 — 30, which becomes sub-dominant at lower z (z ~ 15). 

Going further requires improvements from the observational 
and theoretical point of view. On the theoretical part, the ex- 
istence of PopIII stars with masses of a few tens of solar 
masses is also s uggested b y recent hydrodynamical simulations 
jBrornm & Loebll2003l; iGreif & Bromrnll2006t Ijohnson & Bromml 
' 2006L | 2007l; I Johnson et al.ll2007l ; lYoshida et al.ll2007l ; iGreif etall 
2008 ; Ismifh et alj|2009l) . The explosion of these high mass stars 
disrupts the surrounding environment and delays the formation 
of lower mass stars within the same structure. Merger trees 
dSalvadori et al.l 12007) account for such a mechanical feedback 
which delays PopII/I star formation along a single branch. In a ho- 
mogeneous picture, this delay translates into a reduction of the SFR 
between the two modes. We note that this epoch is very short in our 
model (100-200 Myr). However, the exact evolution of the struc- 
tures are still uncertain. One also needs improved yields provided 
by stellar evolution models at all metallicities and, more specifi- 
cally, SN calculations that are even more uncertain. Note however 
that uncertainties related to the determination of abundances (de- 
rived from ID or 3D model atmospheres) do not change our con- 
clusion on the cosmological importance of PopIII stars. 

Many improvements are possible in the future with regard to 
observations, (i) Additional stellar constraints and abundance mea- 
surements. Better statistics are necessary to construct the MDF at 
very low metallicity, with next generation optical telescopes such 
as GMT. It will then be possible to constrain in a better way the 
epoch of the massive population (whose duration is related the 
critical metallicit y as used in merger trees; ISalvadori et al, I l2007l; 
iTumlinsor] 120061 and the typical mass range of PopIII stars, (ii) 
More complete spectroscopic observations in very metal-poor stars 
are needed in order to obtain a complete set of observations in the 
Strans diagram. Of particular importance is the continued search 
for HMP stars and UMP stars to test their statistical significance 
relative to the bulk of the VMP/EMP stars. (Hi) In future studies, 
one could attempt to better constrain the model using other chemi- 
cal elements, in particular r and s-process elements which can bring 
new constrain ts on the mass range of massive stars, PopIII (e.g. 
lFrebelll2008bl) . (iv) Additional cosmological constraints at high 
redshift. Massive objects may be directly observable, as more and 
more quasars and galax i es are detected at high redshift dFan et all 

120061) and with JWST in the fu- 



l2003l : lKneib et alj|2004l;llve et al 



ture (e.g. IStiavelli & T renti 2009). However, the ex treme brevity 



of the PopIII epoch makes the connection unclear ^Alvarez et al.l 
120071) . A more promising probe would be the observations of high 
redshift gamma-ray bursts, that cor respond, for the longer bursts, to 
the deaths of massive single stars 1 Heg er et al. 2003; Barto rTetail 
120041 ; iBromm & Loebll2006i ; iDaigne et alj|2006f) . Recent observa- 
tions of GRBs at z > 6 (GRB 050904 at z = 6.3, iKawai et al.l 
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6.7, iGreiner et al. 2009) and even 
Tanvir et ajT 20091) are ex- 



8.2, 



l20O6t GRB 080913 at z 
z > 8 (GRB 090423 at 
tremely encouraging, (v) In addition, massive stars would have 
polluted their environment with an initial enrichment of heavy ele- 
ments which could be compared to the one obse rved in the Lyq 
forest along quasar absorption spec tra at z < 6 <SongaUdl200lt 
lAguirre et al.ll 20021; I Aracil et alj|2004h or in the Dam ped Lyman a 
systems ( tLedoux et alj|2003l : IProchaska et al]|2003l) . This requires 
a better understanding of the outflows of metals into the IGM which 
is beyond the scope of this work. In addition, it will be important to 
better understand the role of inhomogeneities on the different pop- 
ulations of the first stars. The uniformity and extent of the metal 
pollution is also under debate, and could in the future be used to dis- 
tinguish between local and re cent pollution and global pollution by 
an earlier population of st ars dMadau et alJl200ll ; lMori et alj|2002l ; 
IWada & Venkatesanll2003l) . 
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